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At the beginning of last century, Gerlach and Lehrer 1 ' 2 observed the rotational Brownian 
motion of a very fine wire immersed in an equilibrium environment, a gas. This simple 
experiment eventually permitted the full development of one of the most important ideas of 
equilibrium statistical mechanics: the very complicated many-particle problem of a large 
number of molecules colliding with the wire, can be represented by two macroscopic 
parameters only, namely viscosity and the temperature. Can this idea, mathematically 
developed in the so-called Langevin model and the fluctuation-dissipation theorem 3 ' 4 , be used 
to describe systems that are far from equilibrium? Here we address the question and 
reproduce the Gerlach and Lehrer experiment in an archetype non-equilibrium system, by 
immersing a sensitive torsion oscillator in a granular system 5 ' 6 of millimetre-size grains, 
fluidized by strong external vibrations. The vibro-fluidized granular medium is a driven 
environment, with continuous injection and dissipation of energy, and the immersed 
oscillator can be seen as analogous to an elastically bound Brownian particle. We show, by 
measuring the noise and the susceptibility, that the experiment can be treated, in first 
approximation, with the same formalism as in the equilibrium case, giving experimental 
access to a "granular viscosity" and an "effective temperature", however anisotropic and 
inhomogeneous, and yielding the surprising result that the vibro-fluidized granular matter 
behaves as a "thermal" bath satisfying a fluctuation-dissipation relation. 

Before describing the experiment in the granular medium we briefly recall the Langevin 
formalism for the Brownian motion of a torsion oscillator immersed in a usual liquid at temperature 
T. This is a canonical problem of statistical mechanics 4 . The Langevin equation of the oscillator is 
I(d 2 &/dt 2 )+a(d&/dt) + G&=C exl (t) + R(t), where d is the angular deflection, /is the moment of 
inertia of the oscillator, a the friction coefficient determining the viscous torque, G the fibre elastic 
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torsion constant, C ext (t) the external torque, and R(t) a randomly fluctuating torque, assumed to be 
a Gaussian white noise of zero mean. This equation is valid on timescales large compared to the 
correlation time of the random torque. 

The Langevin equation can be solved by Fourier transformation. Without the external torque 
( C ext (t) = 0), using the Wiener-Khintchine theorem for stationary processes 4 , one obtains the noise 
Power Spectral Density (PSD), i.e. twice the Fourier transformation of the auto-correlation function 
<6(t)0(t')>, which gives S(0)) = 2q/f 2 ((O 1 o -0) 2 ) 2 + aW\ where 0) = 4g7i is the natural 

pulsation of the oscillator and 2q is the PSD of the random torque R(t) . Moreover, one has the 
relation q = 2ak B T, which guarantees that the system is in thermodynamic equilibrium and that the 
equipartition of energy holds. On the other hand, solving the Langevin equation with the external 
torque C ext (t), and focusing on timescales large compared to the correlation time of the random 
torque, yields the complex susceptibility x( a> ) = ^( co ) /C exl ((0) = Z'(Q))- ix"(a>) with 
X(CO) = l/|/(fi^ - #> 2 )+ iacoY One finds that the ratio S(co)co/{Ax"(co))=k B T is proportional to the 

temperature of the liquid, and is independent of the oscillator characteristics, such as its mass or 
shape. This represents a formulation of the fluctuation-dissipation theorem 4 (FDT), expressed in the 
frequency space. Summing up, measuring the susceptibility and the noise PSD of the Brownian 
motion of the oscillator allows one to consistently determine the viscosity and the temperature of 
the liquid seen as an ideal thermal bath at equilibrium. 

The natural question that arises at this point is how much of this formalism survives if the 
experiment is performed in a strongly non-equilibrium system, notably in a vibro-fluidized granular 
medium 5 ' 6 . A granular system is an assembly of particles, such as sand grains or glass beads, 
interacting by contact forces and featuring a very large number of macroscopic degrees of freedom, 
corresponding in first approximation to the positions and velocities of all grains. However, such a 
system is not in "equilibrium" in the thermodynamic sense, since the thermal energy k B T at room 
temperature is too small to induce any macroscopic grain fluctuations. Nevertheless, due to the very 
large number of degrees of freedom, one expects that an analogy to the thermally-induced Brownian 
motion is possible when the system is externally driven and grain motion occurs by continuous 
injection of energy. 

We have realized the experiment by immersing a torsion oscillator in a granular medium 
composed of glass beads (Fig. 1). The container, filled with the glass beads, is continuously shaken 
by a vertical vibrator, with a high-frequency filtered white noise, cut off below 300 Hz and above 
900 Hz. Notice the important point that a filtered white noise is used in order to obtain a 
homogenous agitation of the granular medium, discarding undesired effects such as pattern 
formation, rolls, and other instabilities in the granular motion 5 . The purpose is not to provide ab 
initio a random torque with white noise spectrum to the oscillator. Indeed, we observe the motion of 
the oscillator in a low-frequency range (10 Hz to 50 Hz) compared to the applied vibration high- 
frequency range. The probe of the torsion oscillator, with various moments of inertia and shapes, is 
immersed at a depth L from the granular surface. The oscillator is otherwise isolated from the 
container and the vibrator, and does not move in the vertical or horizontal directions. An 
accelerometer on the container is used to measure the acceleration spectrum, normalized to the 
acceleration of gravity, A(a>). We quantify the "intensity" of the external shaking by r, defined as 
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the square root of the band power of A(co) in the range about 1 Hz to 10 kHz, i.e. r 2 = J A(f)df 
integrated in the above frequency range, with C0 = 2nf. In the case of a sinusoidal vibration of 
amplitude a s and frequency f s = C0 S /2n , the normalized acceleration spectrum is 

-fj and one recovers the usual definition of the vibration intensity, i.e. 

r= a s O) 2 s / g. For sinusoidal vibrations, T = 1 is the lift-on threshold above which a single grain 
starts to "fly". In this work we employ vibration intensities up to r~ 15. 

We now use the oscillator to measure 1) the noise PSD, and 2) the susceptibility of the 
vibro-fluidized granular medium: 

1) In absence of the external torque ( C ext (t) = 0) the immersed oscillator performs an 
irregular free angular motion, induced by the continuous interactions of the grains against the probe. 
The angular deflections 6 are detected optically using a mirror fixed on the oscillator and, from the 
time-series of the angular deflections, 0(t) we obtain the noise PSD, S((ti), as shown in Fig. 2a for 
several values of r. 

2) While the granular medium is vibrated at a given intensity r, a sinusoidal torque 
C ext ( t) = C e sin( cot) is applied to the oscillator using a permanent magnet fixed on the oscillator and 
two external coils. The complex susceptibility at the given CO is obtained from 
X(co) = 6(co)/C ext ((o) , where 6(co) and C ext (co) are the Fourier transforms of 9(t) and C ext (t) , 
respectively. The amplitude C e of the external torque is small enough to be in the regime of linear 
response, i.e., we measure the linear susceptibility. The complex susceptibility %((o) versus CO is 
obtained by sweeping the frequency of the sinusoidal torque. The modulus of the complex 
susceptibility, \z(co)\, for different intensities of vibration r, is shown in Fig. 2b. 

With both the susceptibility and the noise PSD data, we are now in the position to check 
whether the analogy to the thermal-induced Brownian motion makes sense. First of all, the modulus 
of the complex susceptibility, \z(o))\, is fitted remarkably well by the susceptibility expression for 

-1/2 

, as shown in Fig. 2b. This means that the 



the damped oscillator, \%((o)\ = 



I 2 (col-co 2 ) 2 + a 2 co 2 



deterministic equation of motion, i.e. I{d 2 dl dt 2 ) + a(ddl dt) + Gd=C ext (t) , describes well the 
response of the immersed oscillator, and allows us to define a granular friction coefficient a, or a 
granular viscosity jj.^-a. From the fitting parameters at different r we also deduce that a°^\/ T, 
as shown in the inset of Fig. 2b. 

Secondly, we can obtain the fluctuation-dissipation (FD) ratio S(co)co/ Az"(co) from our 
data. Fig. 3a shows the FD-ratio versus / = CO/27C at different r. The FD-ratio is surprisingly "flat", 
i.e. approximately independent of the frequency in the observed low-frequency range (especially 
compared to what has been measured in other non-equilibrium thermal systems, such as in laponite 7 
and glycerol 8 .) This means that the high-frequency driven agitation of the granular medium acts on 
the oscillator as a source of random torque with white spectrum, at least in the 10 Hz to 50 Hz 
range under investigation. Energy is thus injected at high-frequency, and spreads into a low- 
frequency white spectrum. A roughly flat FD-ratio also provides support for the existence of a 
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fluctuation-dissipation relation in off-equilibrium driven granular steady states. The FD-ratio level 
can thus be used to define an effective temperature, T eff . Fig. 3b shows the averaged FD-ratio levels, 

i.e. k B T eff , versus r. Fitting to a power-law we obtain k B T eff ^r p with p close to 2, thus giving 

Obviously, this is only a first-order approximation. Close inspection of Fig. 3 a apparently 
shows a small frequency dependence. However, an useful torsion signal is only detected around the 
natural frequency of the oscillator, limiting our accessible frequency range. In future work, the 
study of the frequency dependence of the FD-ratio over a large frequency range could give insights 
about possible energy "cascade" in vibro-fluidized granular matter. 

In the frame of this first-order approximation, our results suggest a simply picture. Due to 
the complex dissipation processes between the grains, only a fixed fraction of the energy injected by 
the vibrator is effectively available as granular kinetic energy and is "sensed" by the oscillator. In 
fact, we notice that the order of magnitude of the thermal energy, i.e. k B T eff , as measured here, is 

consistent with realistic values of the mean kinetic energy per particle, as measured by grain- 
tracking methods, e.g., in ref. 9 and 10. Thus, the effective temperature T eff we measure seems to be 

related to the granular temperature, as usually defined in granular gases 5 ' 6 . Notice that the measured 
granular friction coefficient a decreases by increasing the effective temperature, since a °< 1 AT and 
T eff o^r 2 . This is the behaviour of liquids rather than gases. However, the oscillator sees an 

increasing effective temperature and a decreasing granular density as r is increased. A similar 
effect has been reported by Zik et al. (ref. 1 1) by observing the mobility of a sphere immersed in a 
vibro-fluidized granular medium. 

We now study the dependence of T eff on the oscillator properties, as well as on the 
characteristics of the granular medium, such as immersion depth and granular anisotropy. Fig. 3b 
addresses all these issues. The main panel of Fig. 3b shows data obtained using a conical probe with 
different moments of inertia. The upper inset of Fig. 3b shows data obtained using cylindrical 
probes with different sections and diameters. In all these circumstances it turns out that T eff is 

insensitive to the change of these parameters. This is an indication that T e g could be an intrinsic 
property of the granular medium. 

The bottom inset of Fig. 3b shows data measured at different immersion depths, using 
conical and cylindrical probes. Two important aspects needs to be underlined: First, in a layer close 
to the surface the measured FD-ratio shows a strong frequency dependence, and it becomes "flat" 
(as in Fig. 3 a) only at a sufficient immersion depth. (The FD-ratio becomes "flat" about at the depth 
where a minimum appears in the bottom inset of Fig 3b). Underneath the surface layer, we observe 
that T e ff increases with the depth. This is possibly related to a slight decrease of the granular density 

with the depth, as observed in simulations 12 and measurements 13 . Regarding the behaviour of the 
surface layer itself, even though we cannot exclude a relation with density induced features (as 
observed in refs. 10,12 and 13), it ought to be explored more carefully. A second point is the 
difference between T eff measured by conical and cylindrical probes. This is a peculiar granular 

effect and it is probably related to the intrinsic anisotropy of the system under gravity. Although the 
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diameter or the form of the horizontal section of the probe play no role, the angle of inclination of 
the surfaces exposed to the granular agitation, relative to the vertical axis, seems to be a relevant 
parameter. As reported in ref 10, the granular temperature is anisotropic, corresponding to 
horizontal and vertical impulse components. One thus expects in particular the transfer of the 
vertical impulse into torsion momentum to depends on the geometry and the surface state of the 
probe, and this effect will be the object of further investigations. Nevertheless, for a given kind of 
probe, e.g. conical or cylindrical, the present experiment shows that one can satisfy the 
requirements for the Langevin approach to be applicable. 

The most important finding obtained here is that it is possible to separate slow and fast 
degrees of freedom in vibro-fluidized granular matter and to model it, in the very spirit of the 
Langevin approach, in terms of two intrinsic parameters: an effective friction coefficient and an 
effective temperature. In this perspective the oscillator behaves indeed as a Brownian oscillator, and 
the vibro-fluidized granular matter can be seen as a macroscopic "thermal bath". We caution that 
the analogy is only formal, since the granular medium is never in equilibrium in the thermodynamic 
sense. The experiment also shows that the granular matter's mechanical anisotropy and 
inhomogeneity due to the gravitation play an important role, inducing complications that are not 
present in usual gas or liquid media, which can be assumed isotropic and homogenous. 
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LEGEND 

Figure 1 Sketch of the torsion oscillator immersed in a "viscous liquid". Here the "viscous 
liquid" is a granular medium fluidized by strong external vibration (see text for details). The 
granular material consists of glass beads of diameter 1.1 ±0.05 mm contained in a metallic 
bucket of 60 mm height and 94 mm diameter, filled to a height of 31 mm. The oscillator 
has a natural frequency of about 12 Hz, depending on the probe used. The main probe 
used is a cone of apex angle 120° covered with a s ingle layer of glass beads, glued on by 
an epoxy. Probes with different shapes have also been used, in particular cylindrical 
probes with various section and diameters, see Fig. 3. The probe is immersed at a depths 
L . The oscillator allows measurement of both the noise PSD and the complex 
susceptibility. An accelerometer permits to quantify the intensity of the external vibration, 
r. Notice that the oscillator does not move in the vertical or horizontal directions. For this 
reason, below r « 1 the oscillator acts as a fixed-displacement deformation point, and the 
study of the granular dynamics at very low r requires a different approach (ref. 14). 

Figure 2 Noise and susceptibility for different vibration intensities. The probe is conical, 
immersed at about 11 mm from the surface, a, Noise PSD, S(a>), versus the frequency f , 
for different vibration intensities r. The intensity values are, from top to bottom, 11.6, 1 0.0, 
7.3, 5.4, 3.7, 2.2, 1.5 and 1.0. b, Modulus of the complex susceptibility \fl(co)\ versus the 
frequency f = co/27t, for different vibration intensities r. The intensity values are as 
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above. Each curve (continuous line) is fitted (dashed line) with the expression 

-1/2 



\ Z (co)\= l 2 (a? -co 2 j + c?co 2 



with two free parameters, co and a. The moment of 



inertia of the oscillator is / = 1 .5-10 6 kg m 2 , and the applied torque amplitude is 

C exf ('fj = 3.2x10" 5 N m. The parameter a> shows a small dependence on r: it features a 
small increase from 12 Hz to 18 Hz by reducing r. The Inset shows the friction coefficient 
a, obtained from the fit to the curves ^(<o)\ in the main panel, versus r. The dashed line 

is a power law a<*\IT . The l/r dependence is also systematically observed using 
susceptibility data obtained with probes of different shapes. 

Figure 3 FD-ratio and the effective temperature, a, The FD-ratio S(co)(o/(A%"((oj), versus 

f = co/2tt for different vibration intensities r, obtained from Fig. 2. The intensity values 
are, from top to bottom, 11.6, 10.0, 7.3, 5.4, 3.7, 2.2, 1.5 and 1.0. b, The FD-ratio level 
(averaged from 10 Hz to 50 Hz), i.e., the effective temperature T eff , according to the 
expression S((o)(o/iAx"(coj)=k B T eff , versus r, from data in Fig. 3a (black open symbols), 
and from similar data obtained using a conical probe with triple moment of inertia (red 
symbols) and the same L. A power law fitted to the data gives T eff c<r p , with p=2.1. The 

dashed line has equation k B T eff = 3.5x1 0" 10 r 2 . The same FD-ratio levels are measured 
feeding the vibrator by a filtered white noise with a different central frequency (800 Hz, not 
shown). The FD-ratio level depends also on the tribological properties of the granular 
material (especially rich for glass, ref. 15). We have obtained reproductible data using 
aged (i.e., exposed long time to atmosphere) glass beads. Upper inset, the effective 
temperature T eff versus /, obtained using probes of different shapes: a cylinder of circular 
section of 5 mm diameter and polished surface, and a cylinder with "toothed wheel" 
section, with polished surfaces, as sketched. The immersion depth is about 21 mm. The 
same FD-ratio levels are also observed using circular section of 7 mm diameter. Lower 
inset, the effective temperature T eff at r = 9, versus the immersion depth of the oscillator, 
L, for conical (red symbols) and cylindrical (black symbols) probes (see text for details) 
(For clarity, only one symbol every four data values is marked). Notice that the precise 
position of the surface in unknown for a vibrated granular medium, but it is approximately 
at L = (vertical dashed line). The bottom of the container is at L = 31mm. We caution 
that in the layer close to the surface (about 0<L<5 mm), the FD-ratio is anymore "flat", and 
we do not report data in this region. 
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